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ABSTRACT 

Aims. The aim of this work is to determine the multi-thermal characteristics and plasma energetics of an eruptive plasmoid and oc- 
culted flare observed by Solar Dynamics Observatory/ Atmospheric Imaging Assembly (SDO/AIA). 

Methods. We study an event from 03-Nov-2010 (peaking at 12:20UT in GOES soft X-rays) of a coronal mass ejection and occulted 
flare which demonstrates the morphology of a classic erupting flux rope. The high spatial, and time resolution, and six coronal chan- 
nels, of the SDO/AIA images allows the dynamics of the multi-thermal emission during the initial phases of eruption to be studied in 
detail. The Differential Emission Measure (DEM) is calculated, using an optimised version of a regularized inversion method (Hannah 
& Kontar 2012), for each pixel across the six channels at different times, resulting in emission measure maps and movies in a variety 
of temperature ranges. 

Results. We find that the core of the erupting plasmoid is hot (8-11, 1 1-14MK) with a similarly hot filamentary "stem" structure con- 
necting it to the lower atmosphere, which could be interpreted as the current sheet in the flux rope model, though is wider than these 
models suggest. The velocity of the leading edge of the eruption is 597-664 km s~^ in the temperature range >3-4MK and between 
1029-1246 km s~^ for <2-3MK. We estimate the density (in 11-14 MK) of the erupting core and stem during the impulsive phase 
to be about 3 x 10^ cm-^ 6 x 10'^ cm-^ 9 x 10^ cm ^ in the plasmoid core, stem and surrounding envelope of material. This gives 
thermal energy estimates of 5 X 10^^ erg, 1x10^^ erg and 2x10^*^ erg. The kinetic energy for the core and envelope is slightly smaller. 
The thermal energy of the core and current sheet grows during the eruption, suggesting continuous influx of energy presumably via 
reconnection. 

Conclusions. The combination of the optimised regularized inversion method and SDO/AIA data allows the multi-thermal character- 
istics (i.e. velocity, density and thermal energies) of the CME eruption to be determined. 
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1. Introduction 

The solar atmosphere is awash with highly dynamic phenom- 
ena driven by the rapid liberation of magnetically stored energy. 
Magnetic reconnection is thought to facilitate this energy re- 
lease, resulting in solar flares, coronal mass ejections (CMEs) 
and the associated particle acceleration, plasma heating and bulk 
motions. The classic model of solar eruptive events (SEE) is 
the CSHKP model (|Carmichael|1964t|Sturrock|1966l|Hirayama 



1974[[Kopp & Pneuman|1976| ) which con tinues to become more 
ornate over time (for instance |Lin & Eor bes 2000). The basic 
scenario is that of an erupting flux rope (helical magnetic struc- 
ture) which stretches and elongates the coronal magnetic field 
behind it, instigating magnetic reconnection in a vertical current 
sheet. This releases considerable energy, particularly in acceler- 
ated particles that then travel down into the denser lower solar 
atmosphere. Here they lose their energy through collisions with 
the background plasma, producing hard X-ray (HXR) footpoints 
and flare ribbons, heating material that then expands upwards 
filling flaring loops. 

Elares are well observed in the lower solar atmosphere, but 
CMEs have, till recently, mostly been studied at higher altitudes 
(for an overview see [Hudson et al.| ( |2006| ); [Webb & Howard] 
( |2012| )) due to the lack of suitable high cadence observations. 
This means that there is considerable uncertainty as to the spe- 



Send offprint requests 
iain . hannah@glasgow .ac.uk 



cific processes that develop flux ropes and initiate CMEs (e.g 
Demoulin & Aulanier 2010). High in the corona, white-light 
observations typically show CMEs structured as a bright loop 
ahead of a darker cavity containing a bright core. Here the bright 
loop is taken to be due to the increased density of material 
amassed in the leading front of the eruption with the cavity 
likely due to to a flux rope ( [Gibson et al.||2006] ). CMEs have 
been found to evolve in distinct stages, initiation, acceleration 
and propagation (Zhang & Dere 2006) and the initial accelera- 
tion stage of CME synchronize well with the flare HXR ( Kahler| 
et al.|1988t [Sterling & Moore|2005t[Temmer et al.j2008,,2010|) 



and soft X-ray (SXR) emission ( Zhang et al.|2001[ Vrsnak et al.[ 
2005|, while the HXR flux correlates with the acceleration of 
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CMEs ( [Temmer et al.|,201 0). Active regions are significantly 
more likely to be eruptive if they have a sigmoidal morphologly 
( ,Canfield et al.|199 9^ and these struc tures are s een as th e precur- 
sors to flux ropes both observ ational ([Liu et al . 2010; Green et al. 
2011) and in simulations (Torok & Kliem 2003; Ean & Gibson 
,2003, Archontis & Torok 2008 ; Aulanier et al. 2010i r 

Associated with CMEs is the nearby dinmiing of the corona 
during the onset sta ge, seen in soft X-rays ([S terling & Hudson] 
1997]) and EUV (e.g. [Gopalswamy & HanaQka|199 8 ). This could 
be due to the depletion of material as part of the eruption, plasma 
heating or a combination of both. Observations generally point 
to the former: for instance five CMEs studied by H arrison et al] 
(2003 ) using EUV spectroscopic data from CDS showed dim- 
ming due to mass loss and not temperature variations. However 
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03-NOV-10 12:15:02 



Fig. 1. SDO/AIA observations from each of the six coronal filters during the impulsive phase of the 3-Nov-2010 flare when the 
erupting plasmoid is clearly visible in 131 A. The time shown is for the 94A image, with all from 12:15:00 to 12:15:09UT. 



the are also examples of CMEs where the dimming is due to 
plasma heating and not depletion ( Robbrecht & Wang|2010 ). 

Spectral observations of CMEs in the extended corona (> 
I.SRq from Sun centre) with SOHO/UVCS revealed many as- 
pects of CMEs at these altitudes ( |Kohl et aL]2006 ), observing 
spectral lines from CMEs over 0.02 - 6 MK. In one CME a 
"narrow" (lOOMm) structure was observed for 20hrs, with an 
emission measure of about 10^^ cm"^ over 4-6 MK connect- 



ing th e CME core to the flare loops below ( [Ciaravella et al. 
2002| ). This was interpreted as an edge-on view of the current 
sheet and a density of 5 x 10^ cm"^ was estimated using a 40Mm 
line-of-sight component. Similarly hot current sheets have been 



observed in another long-duration event (Ko et al. 



well as faster ones ( [Lin et al.||2005| ) where the leading edge 



2003) as 



and CME core were moving at 1 939 and 1 484 kms" respec 
tively. Hotter current sheets have been observed in some CMEs 
as well (Innes et al.^ 20()3{ [Bemporad et al.||2006^ Schettino] 
|er^|2010| ). |Landi et al.| ( |2Q10| ) used Hinode^IS to determine 
the Diff'erential Emission Measure (DEM, see @ of the CME 
lower in the solar atmosphere (at IARq). They found a multi- 
ple peaked DEM, with two components at 0.13 MK of density 
2Q10-11 another at 0.5 MK of density 10^-^ cm"^ trav- 

elling together. The emission at about IMK was due to the line- 
of-sight/background corona however SXR emission was also de- 
tected by Hinode/XRT, enveloping the core and suggestive of 
material at 5-10 MK. The analysis of the post eruption arcade 
found similar plas ma to the current s heet, both coronal plasma 
heated up to 3 MK ( |Landi et al.|2012| ) 

High spatial (~ 1.5'') and temporal resolution (~ 12 s) EUV 
full disk images in multiple wavelengths is now available from 
the Atmospheric Imaging Assembly (Lemen et al. |2012| ) on 
the Solar Dynamics Observatory (SDO/AIA), with several au- 
thors taking advantage of this data to study CME formation and 
propagation from the low corona (e.g. |Patsourakos et al.||2010[ 
Kozarev et al. 2011| ). One particular event, which peaked in 



GOES SXRs at 12:20UT 03-Nov-2010 (Figure[T]), of an occulted 
flare and erupting plasmoid, shows unambiguously the forma- 
tion of a flux rope and has subsequently been studied in detail 
dReeves & Golub 201 1! [C heng et al.|201H[Fbullon et al.|201H 



Bain et al. 2012; Savage et al. "201 2; 'Cheng et al.'2012), and wiU 



be studied in detail in this paper. This event shows the formation 
of a flux rope clearly in 131 A, images with a peak in their sen- 
sitivity at IIMK, though dimming at the core in channels sen- 
sitive to cooler emission (171 A, 211 A), shown in Figure [l] As 
the plasmoid move outwards, a filamentary, stem-like structure 
forms behind it, again sharply defined in 131 A, suggestive of a 
large-scale current sheet (Reeves & Golub 2011 ; Cheng et al.[ 
2011 [[Savage et al.|2012| ). The kinematics of the eruption in the 



diff'erent channels has been found to be slower (about 600 kms ^ 
versus 1 400 kms~^) in those sensitive to hotter emission (Cheng 
[et al.|2011l[Bain et al. 12012"] ). 

The propagation of a metric Type II radio burst ahead of 
the CME in EUV is faster still (nearing 2 000 kms"^) suggest- 
ing a piston-drive shock, in which it is moving too quickly 
for the ambient plasma to flow behind the driver ( [Bain et aT 
|2012). This work also found that HXR emission observed with 
RHESSI had two non-thermal (18-40 keV) sources, one compact 
low in the corona, the other extended and high in the corona, 
close to the erupting plasmoid. The later could be due to elec- 
trons accelerated in the stem/current sheet. Although this event 
superficially matches the cartoon model of flux rope eruption 
things are more complicated: jFoullon et al. ( [201 1) found Kelvin- 
Helmholtz waves along the top edge of the erupting plas ma, pos- 
sibly allowing a secondary reconnection process. Savage et al. 
( [2012| ) studied the inflow and outflow into the CS region, find- 
ing the hottest material inflowing fastest and correlating with 
the RHESSI thermal X-ray emission. Pairs of inflowing and out- 
flowing material was found to be travelling away from the region 
faster than they arrived, suggesting acceleration. 
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Since SDO/AIA provides six different wavelength channels, 
sensitive to a range of coronal temperatures, they can be used 
to determine the DEM of the different parts of the CME. Cheng 
et al. ( |2012| ) did this for this event, as well as two others, produc- 
ing DEMs for different CME regions: the core of the erupting 
flux rope and the leading front. They found the cores to be have 
emission over a broad distribution of temperatures (3 < T < 
20MK), with an average temperature of about > 8MK and den- 
sities of about lO^cm"^. The leading front of the eruptions were 
cooler (about 2MK on average) with a narrower temperature dis- 
tribution (1 < r < 3MK) similar to the pre-eruption coronal 
material but with slightly higher densities. They interpreted the 
leading fronts to be the showing signatures of compressed ambi- 
ent coronal plasma. Dimming, about 10 minutes after the erup- 
tion, was found to be due to the rarefaction of material. 

In this paper we use a different approach, regularized in- 
version, to produce the DEM not for a few regions within the 
CME but for every single pixel in the set of six SDO/AIA im- 
ages, every 12 seconds throughout the event. Thus allowing the 
highly dynamic multi-thermal properties of the CME to be inves- 
tigated. The regularized approach detailed in this paper (seeQ 
is an optimized version of our original code, introduced, tested, 
and compared to th e approach used by |Cheng et al.| ( |2012| ) in 
Hannah & Kontar| (|2012 ). The version present here is able to 
compute 10 000s of DEMs per second, calculating both horizon- 
tal and vertical and uncertaintie s (se e |2.2| ) and without forcing 
a prescribed DEM model. In ^3.2 



the properties of the DEM 
in different pixels of the event are d iscussed, with the maps and 
their uncertainties presented in ^33 We then investigate the time 
evolution of the DEMs in terms of the morphology, total emis- 
sion and propagation speed for the different temperature ranges, 
in ^3.4| As the erupting plasmoid and filamentary stem struc- 
ture is well defined at high temperature (11-14MK), we use the 
emission measures in ^33] to estimate the densities and thermal 
energy of the different components. 

2. Differential Emission Measures 

2.1. Introduction 

The signal detected g(x,yjj) in each pixel (x,y) in the i^^ 
SDO/AIA filter at time t can be interpreted as 



X 



K(T, i)^(T, X, y, t)dT + Sg(x, y, t, /), 



(1) 



where K(T,i) is the temperature response of the f'^ filter 



(Q'Dwyer et al. 2010 Boerner et al. 2012). The Differential 



Emission Measure DEM ^(T, x,y, t), at a particular instance in 
one pixel, is taken to be ^(T) = n(h)^dh/dT [cm"^K"^], which 
is the electron density along the line-of-sight h of the emit- 
ting material at temperature T. The inherent noise in the obser- 
vations Sg(x,yjj), due to counting statistics, background and 
instrumental uncertainties, will be significantly amplified when 
a direct solution of Eq. ([T]) is attempt ed. This recovery of the 
DEM is an ill-posed inverse proble m ([Tikhonov||1963[ |Bertero 
eraL][T985| jCr aig & Br own 1986^ |Prato et al.||2006] ) with the 
growth of uncertainties resulting in spurious solutions. To re- 
cover useful solutions of the problem given in Eq. ([T]), it needs 
to be constrained to avoid the noise amplification in the DEM 
solution. The constraints are often physically motivated addi- 
tional information a priori known about the solution. An ad- 
ditional challenge here is that we need an approach that can 
quickly and robustly calculate DEMs since there is a wealth of 



SDO/AIA data to study this highly dynamic eruption. The sim- 
ple isothermal/ratio approach (e.g. Vaia na et al.||1973| ) places a 
very strict, though possibly erroneous, constraint on the DEM 
solution and so is ill-suited to this problem. An alternative is to 
forward fit a chosen DEM model to the data, for example with 
mult i-Gaussians (e.g. |Aschwanden & Boerner||2011| ) or splines 
(e.g. [Weber et al.||2 004), the latter used for this CME event by 



Cheng et al. (2012|. Again, if the assumption about the form of 



the DEM is wrong then the results are difficult to reliably in- 
terpret. These techniques can be computationally slow, only al- 
lowing select regions to be analysed, not taking full advantage of 
the SDO/AIA data. Other approaches, such as MCMC ( Ka shyap] 
& Drake 1998 ), can produce very robust results, without requir- 
ing predefined functional forms to the DEM solution, and uncer- 
tainty estimates but are computationally slow and so are again 
inappropriate to obtain DEM maps. 

Tikhonov regularized inversion recovers the solution by lim- 
iting the amplification of the uncertainties through assumptions 



about its ''smoothness" by using line ar constraints ( [Tikhonov 
1963} jCraigjpWTl [Prato et al.|[2006l ). This approach does not 



make any assumption on the functional form of DEM and we 
have implemented an algorithm ( Ha nnah & Kontar|2012[ ) to re- 
cover the DEM from solar data using Tikhonov regularization 
and Generalized Singular Value Decomposition GSVD ( [Hansenj 
1992 ). This code, originally developed and extensively tested for 
the inversion of X-ray data ( Kontar et al.|200 4, 2005 ), naturally 
provides estimates to the uncertainties in the solution, an im- 
portant feature, when trying to recover a large number of DEMs. 
Recently, a simpler and computationally fast (computing 100s to 
1 ,000s DEMs per second) S VD regularisation approach has been 
applied to AIA data (Pl owman et al.^201 2 ), however it does not 
produce error estimates. Regularization has also been used with 
solar data to determine the 3D emission measure through the use 



of solar rotation or multi- spacecraft observations (Frazin et al. 



|2005 , 2009 ). In general, regularized approaches do not require 
a positive solution but one can be obtained by using an iterative 
method such a Reiterated Conjugate Gradients, which was used 
for solar data by Fludra & Sylwester] ( |1986j ). Our approach uses 
one iteration to help with the weighting, and positivity, of the 
final solution whilst maintaining the fast computational speed. 

2.2. Regularization Method for DEM maps 



The full details of the method are in [Hannah & Kontar[ ( [2012| ) 
and this previous version of the algorithm allows a variety of dif- 
ferent controls to the nature of the regularization but is too com- 
putationally slow to practically produce DEM maps and movies 
(running at about a few DEMs per second). We have produced an 
optimized and simplified version of the softwar^ detailed be- 
low with considerable speed gains. As the DEM calculation per 
pixel location is independent of other pixels in the same image, 
further speed gains were achieved by parallization of the code 
within IDL using Bridges. The results is that we can achieve 
about 1 000s DEMs per second on a standard multi-core cpu lap- 
top/desktop and 10000s on multi-cpu machines. 

To obtain the DEM maps we perform regularization with the 
zeroth-order constraint and no initial guess solution. For the de- 
tailed explanation of these para meters and the affect they have 
on the regularized solution see |Hannah & Kontar[ p012| ). The 
specific problem to be solved here is 



K^-| +^||L^|p = min, 



(2) 



Available online: http://www.astro.gla.ac.uk/~iain/denireg/map/ 
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where the constraint matrix L is the identity matrix (normalised 
by the temperature bin width) and the tilde indicates normalisa- 
tion by the error in the data Sg. The GSVD singular values (y/, 
Pi) and vectors (u, v, w) produced from the GSVD of K and L 
provide a ready DEM solution as a function of the regularisa- 
tion parameter A ( Hansen 1992) , i.e. 



0? /(g-U/)W/ 



7i 



(3) 



where 0/ = yi/fii. To find the regularization parameter A that 
produces a solution that matches the required - ^ data 
space) we use coarser sampling of A than in the original soft- 
ware which helps with the computational speed increase. Even 
with this approach we achieve of typically between 0.5 and 
1.75 with our regularized solutions, which are identical (within 
the error bars) to the results from the previous code. The above 
procedure is run twice (as in the original code) with the result 
from the first run ^x(T) being used to weight the constraint ma- 
trix as Y?.. = dlogTj/^;i(Tj) before the second GSVD is per- 
formed and final solution found. The horizontal errors are 
again taken to be the spread of the resolution matrix R^K for 
each temperature bin ( [Hannah & Kontar|2012| ), where is the 
regularized inverse assumed to be close to the generalized in- 
verse i.e. ^ R^g. For the vertical error w e do not use the 
Monte Carlo approach as in [Hannah & KontS1(2Q12) but take 



it to be A^iTY =< S^(Ty >= ||Ri^g||^ This produces a similar 
result but is considerably quicker to compute. 

3. DEM Maps of 03-Nov-2010 Eruptive Event 

3.1. Data Reduction 

The data between 12:12:02 UT and 12:19:02 UT (with a set 
of six diff'erent channel images every 12 seconds) were pre- 
pared to level 1.6 by deconvolving the point spread function us- 
ing aia_deconvolve_richardsonlucy . pro before being pro- 
cessed with aia_prep.pro. The co-alignment of the images 
were checked using aia_coalign_test .pro and found to be 
sub-pixel accurate. The resulting images were reduced down to 
sub-maps covering -1150''to -850''in the x direction and and - 
500''to -200''in the y direction, 501 by 501 pixels. 

For each pixel position (x,y) across the six diff'erent fil- 
ters (/ = 1 ... 6) in the 12 second set t provides the signal 
g(x,yj,i), in units of DN s"^px"^ from which the associated 
errors Sg(x, y, t, i) are calculated using the readout noise and pho- 
ton counting statistics. The set of six values, and their associated 
errors, with the temperature response functions are passed to the 
regularized DEM map code, as detailed in ^ The SDO/AIA 
temperature response functions are calculated using both the 
CHIANTI fix (an empirical correction for missing emission 
in the CHIANTI v7.0 database in channels 94A and 131 A) 
and EVE normalisation (to match EVE spectroscopic full-disk 
observations), further details available online |^ We use fairly 
broad temperature binning, twelve pseudo-logarithmic bins be- 
tween 0.5 to 32MK, as to avoid a large number of tempera- 
ture bins slowing down the code. These are larg er tha n the de- 
signed achievable temperature resolution ( Judge 2010| ), yet nar- 
row enough to recover test models , from testing the recovery of 
Gaussian model DEMs, as in Hannah & KontarX2012j . 



SSW/sdo/aia/response/chiantifix_notes.txt 



3.2. DEM from individual pixels 

The resulting DEMs found for eight distinct pixel locations in 
the images from Figure [T] are shown in Figure [2] These pixels 
are shown for the 131 A and 211 A images but the data from all 
six coronal-filter images were used to recover the DEMs. The 
vertical and horizontal errors shown are calculated as detailed in 
^2.2| The vertical error is just the linear propagation of the uncer- 
tainty on the source data but the horizontal error is more com- 
plicated as it gives a measure of the temperature spread (from 
the resolution matrix) as well as the quality of the reg ularized 
solution (i . e. stro ng off'-diagonal terms) in each bin, see Hannah 
|& Kontar] ( |20T2l ) for further details. The DEMs of the different 
positions show a variety of forms but clearly none are isother- 
mal and the majority cannot be represented by a single Gaussian 
model. Shown for comparison are the DEMs obtained in these 
pixels before the eruption at 12:12:02 UT. 
For each specific pixel we have: 

Pixel 1 - Core of the erupting plasmoid: There are clearly 
two distinct components, one peaking in the range 1.5-3MK, 
the other at 8-14MK. The lower temperature one is likely 
the background corona along the line-of- sight whereas the 
hotter from the plasmoid itself. The pre-eruption DEM has 
no high temperature component; 

Pixel 2 - Filament/stem behind the plasmoid: Again two 
components at similar temperature ranges to the core of 
the plasmoid but the low temperature component has larger 
emission (due to higher densities lower in the corona 
compared to Pixel 1) and the higher temperature one has 
considerably weaker emission (due to a smaller amount 
of heated material accumulating in this location compared 
to the core. Pixel 1). The pre-eruption DEM has no high 
temperature component; 

Pixel 3 - High corona away from the event: A single compo- 
nent covering 1-3MK with expected weak emission, since 
lower density higher in the corona, and identical (within the 
errors bars) to the pre-eruption DEM; 

Pixel 4 - Corona away from the event: A broader DEM with 
slightly higher peak emission than Pixel 3, likely due to it 
being lower in the corona but again peaking over 1-3MK, 
and similar (within the errors bars) to the pre-eruption DEM; 

Pixel 5 - Corona near the event: Higher emission and broader 
DEM than similar location of Pixel 4 indicating the presence 
of more hot coronal loops. The emission at high temperatures 
(above 6MK) is zero within the errors bars. The pre-eruption 
DEM is slightly higher, though consistent within the errors, 
so may indicate a slightly hotter and denser loop(s) moving 
out of this region; 

Pixel 6 - Low corona flare emission: Very high emission 
across all temperatures with a bright peak at 1.5-3MK 
and 8-14MK for the flare heated material expanding up 
above the limb to be visible. The pre-eruption DEM shows 
emission across 6-8MK and 8-1 IMK but weaker at higher 
temperatures, which might be due to the behind-the-limb 
flare; 

Pixel 7 - Envelope just ahead of the plasmoid: Weaker emis- 
sion than the core of the plasmoid (Pixel 1) but again a sim- 
ilar DEM with two components, peaking at 1.5-3MK and 
8-14MK, the latter likely due to the accumulation of hotter 
coronal emission in this leading front of the eruption. The 
pre-eruption DEM does not show this hotter (8-14MK) emis- 
sion; 

Pixel 8 - Further ahead of the plasmoid: Similar to the low 
temperature component of Pixel 7 (which is physically closer 
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Fig. 2. (Left) The 131 A and 21 1 A images from Figure[l]with eight pixels highlighted. (Right) The corresponding regularized DEM 
solutions ^(T) (using twelve temperature bins) and the vertical A^(T) and horizontal A log T errors (the maximum of the error from 
the resolution matrix or the temperature bin width) for the highlighted pixels using the data from all six image in Figure [T] The 
dashed lines show the DEM solutions in those pixels before the eruption at 12:12:02 UT whereas the solid lines are from 12:15:02 
UT. 



to the core of the erupting material) but no higher tempera- 
ture component. Consistent (within the error bars) to the pre- 
emption DEM. 

Overall the DEMs in Figur e [2l present a simila r scenario to that 
found by previous authors ( [Cheng et al.|2011[ |Bain et al.|2012| ) 
who were studying the behaviour just in the different SDO/AIA 
channels, as well as the DEM analysis from Cheng et al. (2012). 
Given the inherent difficulties in recovering DEMs this is reas- 
suring that their bulk features are consistent with the original 
SDO/AIA data. Here we have a hot (8-11, 11-14MK) plasmoid 
core, well-defined with a filamentary stem-like structure behind 
it, similarly hot. This connects down to the flare loops filling 
up with even hotter material (14-19MK), which have expanded 
just above the limb, consistent with the classic flux rope eruption 
model. At cooler temperatures there is strong emission ahead of 
the erupting plasmoid but there is also a line-of- sight coronal 
component (peaking about 1.5-2MK) in all the pixels chosen. 

3.3. DEM Maps 

To get a better understanding of the spatial structure of the event 
we can create DEM maps for each temperature range, shown 



in Figure [3] Many pixel locations in these maps show emission 
across a broad range of temperatures, highlighting the complex, 
multi-thermal nature of the line-of- sight DEMs. The tempera- 
ture maps for 0.5- IMK and 8-1 IMK do resemble the images 
from the 171 A and 131 A channels respectively showing that the 
emission here is mostly from one component of the channels' re- 
sponse. However for the other temperature maps, features start to 
appear that are not apparent in the SDO/AIA images themselves, 
such as there being little 4-6MK emission about the core of the 
erupting plasmoid. This emphasizes the complicated tempera- 
ture response of SDO/AIA filters and the necessity of the DEM 
inference to interpret the temperature or density structures in the 
corona. Before discussing these spatial properties it is prudent to 
investigate which of the regularized DEMs are robust solutions. 

One major advantage of our regularization approach is that 
it produces information about the vertical A^(T) and horizon- 
tal A log T uncertainty in the regularized solution (see |2.2| ) and 
so we can plot uncertainty maps, shown in Figure |4] The ver- 
tical error A^(T) (top row, given as a percentage error of the 
DEM signal) show the largest uncertainties are predominantly 
from the regions where the emission measure is very low as 
well as at the highest temperatures. The latter is due to a com- 



5 



Hannah & Kontar: Multi-thermal dynamics and energetics of a coronal mass ejection 







1111111^ ^ ' 'fl 


1.5-2 MK 

1 








^^^^ 4^ ' ^^^1 










^^^^^^H 19.5 


20.5 2 


19-25 MK 

i * 


25-32 MK 

\ 

> 



Fig. 3. DEM maps in different temperature ranges for the 3-Nov-2010 flare at 12:15:02 UT. The same pixel position in the twelve 
maps corresponds to the regularized DEM solution ^(7) (as shown in detail for eight pixel positions in Figure [2]) to the set of data 
values from the corresponding pixel in the six images shown in Figure [l] The colour scale is log^(r) and is in units of cm"^K"^ 



bination of SDO/AIA having lower sensitivity to the temper- 
atures > 20MK and there also being little plasma emitting in 
this temperature range. Hence the data, and also DEM, are very 
noisy. Conversely, at the temperatures where SDO/AIA has a 
very strong response (about 0.5 to 3MK) and there is bright solar 
emission, the signal, and hence DEM, has a very small associ- 
ated vertical error. 

The horizontal error in the solution (bottom row Figure |4]) 
are mostly slightly larger than the SDO/AIA minimum theoret- 
ical limit of temperature resolution A log T = 0.1 ( Judge 2010 ). 
The largest uncertainties again come from regions where there 
is a weak signal and/or noisy data, this is particularly evident for 
the largest temperatures. But there are are also positions where 
the data is noisy in the hottest temperature bands yet the hori- 
zontal error is small. This is due to the regularization approach 
robustly determining that there is very weak/no emission at these 
positions and temperature. Where the horizontal errors (temper- 
ature resolution) are small it means that the emission is confined 
to that temperature range but with large errors (> 0.3) it suggests 
the emission shown in that temperature band could actually be 
attributed to a neighbouring temperature bins. Therefore at tem- 
peratures where the horizontal uncertainty are large the DEM 
solution should be treated with caution (see [Hannah & Kontar| 
2012) and additional analysis of the residuals could be useful 
(e.g. |Piana et al.|2003| ). 

Using this error information we can restrict the DEM maps 
to only show the values with the smallest vertical and horizontal 
errors and this is shown for A^(r) < 30% and A log T < 0.25 
in Figure [5] At the hottest temperature (25-32MK) there is vir- 
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rived from the SDO/AIA observations for the chosen tempera- 
ture ranges (colored lines). 



tually no emission due to the largest uncertainties and poor con- 
fidence in the DEM solution and any weak emission shown is 
likely unreliable. In 14-19MK the emission is predominantly 
coming from the flare heated material that has risen above the 
limb though there is a hint of very faint emission from the plas- 
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Fig. 4. Maps of the (top row) vertical uncertainty A^(r) and (bottom row) horizontal uncertainty, or temperature resolution A log T, 
in the regularized DEM solution shown in Figure [3] for four temperature ranges. 




Fig. 5. DEM maps in different temperature ranges for the 3-Nov-2010 flare at 12: 15:02 as shown in Figure |3|but only showing pixels 
where A^(r) < 30% and A log T < 0.25. The colour scale is log ^(7) and is in units of cm"^K"^ 



moid core. These loops show the strongest emission in 8-1 1 and 
11-14MK but also in these temperature ranges is the plasmoid 
core, an envelope of surrounding material and the filamentary 
stem behind it, all sharply visible. The plasmoid core has slightly 
stronger emission in 11-14MK than 8-1 IMK. Only a single thin 
stem structure is clearly defined in 11-14MK but two, in a V- 
shape diverging away from the plasmoid core, are visible in 8- 
1 IMK. In both cases their high temperatures, and the inflow re- 



gion to these stems, is not explained by the classical flux rope 
eruption model. Some additional heating maybe present and the 
Kelvin-Helmholtz waves to the top of the erupting material are 
a signature of a reconnection layer surrounding the erupting flux 
rope, producing a secondary means of releasing energy ( ,Foullon| 
[itaL 2011, 2012). 

At cooler temperatures (<6MK) the core of the plasmoid dis- 
appears, which is likely the cause of the dimming observed in 
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Fig. 6. Emission Measure maps for various times in 2-3MK, 6-8MK and 8-1 IMK (top, middle, bottom rows). These maps have 
been rotated by 110°. The second column is of the time interval shown in previous Figures. A movie of this is available with the 



electronic version of this paper, covering the 1-1.5 to 14-19MK from 12:12:02 UT to 12:19:52 UT, see Figure A.l 
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Fig. 8. Running differences of the DEM maps (^(7, jc, - ^{T,x,yJi-i)) for 2-3MK (top row) and 11-14MK (bottom row). 
The leading edge of the erupting plasma was tracked and is indicated by the red horizontal lines The resulting height profiles as a 
function of time for these and other temperature ranges is shown in Figure |9] and 10 These maps have been rotated by 110°. 



some SDO/AIA images but it is not clear whether this is due to 
heating or mass loss. To try and understand this we look at the 
time evolution of the maps in |3.4[ The V-shaped structure of 
the stems is still present though, clearly in 6-8MK but slightly 
obscured by the bright surrounding coronal emission in 3-4 and 
4-6MK. Ahead of the plasmoid core we start to see the lead- 



ing edge of a front of material, visible in 2-3, 3-4 and 4-6MK. 
This density increase of material is likely due to a compression 
front being driven ahead of the erupting plasmoid. At the coolest 
temperatures (0.5-1, 1-1.5 and 1.5-2MK) there is some faint ev- 
idence of the leading edge moving out but it is dominated by 
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Fig. 9. Height (above an arbitrary origin) versus time profiles of 
the leading edge of the erupting plasma for each of the tempera- 
ture ranges the propagation could be reliably tracked. 



the emission from several hot coronal loops associated with the 
occulted active region and flare. 

3.4. 7//T?e Evolution 

As the 03-Nov-2010 solar eruption is highly dynamic we investi- 
gate the DEM maps as a function of time, which we have shown 
in three temperature ranges at 1 minute intervals in Figure [6j A 
movie of this is available with the electronic version of this pa- 
per, covering 1-1.5 to 14-19MK from 12:12:02 UT to 12:19:52 
UT, see Figure [ATT] We have also converted the DEM maps into 
EM maps by removing the factor In lOrdlog T, resulting in the 
values shown being EM = n^h [cm"^], making it easier to calcu- 
late densities (see ^3.5| ). Again we have only shown the emission 
in the pixels and temperature ranges that have A^(T) < 30% and 
A log T < 0.25, as in Figure [5] 

In the 2-3MK temperature range (top row Figure |6]) we see 
the leading edge of the eruption clearly at 12:15:02 UT but af- 
terwards gone and the loops bending/collapsing back. There is 
some strong emission ahead of erupting flux rope at 12:14:02 
UT but that disappears at later times, likely erupted. Lower in 
the corona the loops barely change over time and only the high- 
est ones modify, being bent back after the eruption, leaving a 
void where the stem forms in higher temperatures. In 6-8MK 
(middle row Figure [6]) we see the initial eruption of the plasmoid 
core and stems behind it at 12:14:02 UT but at later times there 
is a void of material at this temperature, again suggestive of ma- 
terial loss (or spreading out) due to the eruption mass loss. There 
is increased emission from the loops that have bundled up below 
the eruption to the right indicating possible heating of material 
in these contracting loops. 

The erupting plasmoid and the filamentary stem structure 
connecting it to the lower solar atmosphere is sharply defined in 
8-1 IMK (bottom row Figure [6]). The strong emission from the 
core of the plasmoid quickly weakens as it moves higher leaving 
the strongest emission from the vertical stem behind it. It would 
be tempting to immediately classify this structure as a current 
sheet, given the superficial similarity to the classical eruption 
model. However this stem is relatively thick and appears to be 
made up of several thinner structures. This maybe a projection 
eff'ect with multiple thin stem misaligned along the line-of-sight. 



A bigger issue is that this structure shows strong emission in this 
hot (11-14MK) temperature range and again the classical model 
is unable to explain the presence of such hot material inflow- 
ing to the current sheet. Below to the right of the eruption (or 
to the top right in the original orientation) hot loops bundle up 
between 12:16:02 UT to 12:18:02 UT This is associated with 
Kelvin-Helmholtz waves (Foull on et al.|20lT]) and p ossible den- 
sity increases due to outflowing material (Savage et al. 2012). 
In the classic eruptive model this region would be associated 
with shrinking loops and the outflow from the reconnection re- 
gion. However the dynamics here are clearly more complicated 
and difficult to unambiguously ascertain due to the strong emis- 
sion lower in the atmosphere. In the 12:17:02 UT and 12:18:02 
UT images another large coronal loop develops to the left of the 
main eruption. This fainter emission in 11-14MK is likely flare 
heated material that has expanded into this large loop. STEREO- 
B/SECCHI images from behind the limb do show a separate core 
flaring region located at the base of this hot loop ( Savage et al.| 
[2QT21 ). 

By summing over the whole of the region shown in Figure |6] 
we can produce a time profile in each temperature band, shown 
in Figure 7]in comparison to the GOES SXR emission. We have 
computed this between 12:12:02 UT and 12:19:50UT (the time 
range shown in the movie available in the electronic version of 
this paper, see Figure |A.1[ ) as saturation from the flaring loops 
starts to dominate and produces unreliable EM maps in these 
regions. From Figure [7] we see that there is a slight increase in 
most of the temperature ranges as the eruption begins (12: 13UT) 
which then drops off as the material is ejected (12:15UT), sug- 
gesting this initial change is due to the heating, then the loss of 
coronal plasma. After this the hottest temperature ranges (6-8, 8- 
11 and 11-14MK) increase in a similar manner (though slightly 
delayed) to the SXR emission observed by GOES. This is ini- 
tially due to further heating in the current sheet and the shrink- 
age of loops below the eruption with heated material outffowing 
the reconnection region. The later sharp increase is likely due to 
the bright post-ffare loops rising above the limb. At the cooler 
temperatures it will be later in the event before the current sheet 
and ffare heated material (if visible above the limb) cool into 
these temperature range. This explains the decrease in the emis- 
sion at these temperatures after the loss of material during the 
eruption. The very small decrease at the smallest temperatures 
(0.5-1, 1-1.5, 1.5-2MK) suggest that this emission is dominated 
by the background and line-of-sight coronal emission over the 
whole region. 

To investigate the time profile of the erupting material in dif- 
ferent temperature ranges we have manually tracked the front of 
the emission in running difference maps of the emission mea- 
sure, i.e. ^(r, x,y, ti) - ^(r, X, J, ti-i). This is shown in Figure [8] 
for 2-3 and 11-14MK, with black indicating a decrease in emis- 
sion, white an increase and the colour lines showing the manu- 
ally tracked front of the eruption. This type of analysis has been 
performed before but in individual SDO/AIA wavelength chan- 
nels ( [Cheng et al.||201l[ Bain et al. |2012 ) and not for specific 
temperature ranges. The former has the advantage of a clearer 
time evolution as the channels are a ^3s duration image every 
12s whereas the EM maps are effectively 12s cadence due to 
the use of six different channels in their calculation. This makes 
the determination of the leading edge tricky and the resulting 
height (above an arbitrary position) as a function of time for 
the different temperature ranges are shown in Figure |9] Here we 
have linearly fitted the data to produce a velocity as a function 
of temperature. These velocities and there associated error (the 
largest between the error in the fitted parameter and from three 
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Fig. 12. (Top) Emission Measure maps for various times in 11-14 MK. (Bottom) The emission measure profile at various heights 
for one minute time intervals in 11-14 MK, indicated by the dotted lines in the top panel. 
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Fig. 10. Velocity of the leading edge of the eruption as a function 
of temperature, obtained from linear fits to the data shown in 
Figure |9] The vertical uncertainty is taken as the largest of the 
error in the parameter from the linear fits shown in Figure |9] and 
three AIA pixels travelled in 12s. 




Fig. 11. Emission Measure map for 11-14MK at 12:15:02 UT 
with the diff'erent components identified. 



SDO/AIA pixels travelled in 12s) are shown in Figure 10 In the 
cooler temperature range the leading edge of the eruption travels 
at 1 246 kms"^ in 0.5-lMK decreasing to 1 029 kms"^ in 2-3MK. 
This general trend of a slower eruption at higher temperature is 
more evident above 3-4MK once we are in the range associated 
with the plasmoid eruption itself. From 3-4Mk to 11-14MK the 
velocities are consistently between 597 kms"^ to 664 kms"^ sug- 
gesting the multi-thermal structure is moving outwards together. 
These two distinct structures at diff'erent velocities are consistent 
with those found by tracking the eruption i n the SDO/AIA chan- 
nels (|Cheng et al.||201l| [Bain et al.||2012| ), with the 131 A and 
335A at 667 and 747 kms"^ and the cooler 211 A and 193A at 



1 154 and 1 439 kms'^ ( |Bain et al.|20T2l ). 



3.5. Density and Thermal Energy 

To recover the density from the EM maps an estimate of the line- 
of-sight depth h of the emission is required. As the structures are 
sharply defined at higher temperatures (visbile in Figure [TT] la- 
belling the three distinct components: plasmoid core, envelope 
and stem) we can use their apparent width in the EM maps as 
the depth assuming spherical/cylindrical symmetry. We focus 
just on 11-14MK here as the structures are similar in shape to 
those in 8-1 IMK (bottom row of Figure |6]) and at lower tem- 
peratures this would be difficult as the shape of the structures is 
far harder to estimate reliably. The density of the sheath of the 
Kelvin-Helmholtz waves has been estimated at lower tempera- 
tures (about 4MK) using emission measures found using both 
our regularized approach and Gaussian model forward-fitting 
techniques, findin g similar values, an d is the subject of a sep- 
arate publication ( [Foullon et al.|2012} . 
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Table 1» Thermal parameters for the different features shown in Figures] 11 |andp^ 
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The EM maps in 11-14MK for 5 time intervals during the 
plasmoid eruption are shown in Figure [12] (top row). Here we 
have plotted the EM as a function of position (bottom row) for 
various slices through the erupting structure, indicated in the EM 
maps by dotted horizontal lines. From these slices we can deter- 
mine a width of the structures of interest: the erupting plasmoid 
core, envelope and stem. Assuming a spherical/cylindrical ge- 
ometry in these structure we can use the width as the depth, giv- 
ing the line-of-sight extent h. Since EM - n^h this readily gives 
an estimate of the density in the structure and these values for the 
various structures at the different times are given in Table[T] Here 
we see that the density in the core region (about 3 x l(rcm"^) 
does not change between the two time intervals shown before 
it fully erupts. The envelope of material surrounding the plas- 
moid core has a lower density (3 x lO^cm"^) but still shows a 
sharp transition to the background corona as that has effectively 
no emission in this temperature range. 

The density in the stem region does decrease over time 
but that is initially due to a widening of the structure (seen at 
12:16:02 UT onwards) rather than a reduction in the emission 
measure. As discussed before, it is difficult to attribute such a 
wide hot structure to the monolithic current sheet as proposed 
by the classic eruption model. If we instead assuming the struc- 
tures are more elongated, using a line-of-sight component of 
h = lOOMm (shown in Table [T]), we find that the density of the 
stem does not vary much with time. In fact these slightly smaller 
densities (about lO^cm"^) are fairly similar for both the core and 
stem regions. 

The thermal energy of this emission can be calculated via 
Wt = 3nVkBT, where we take T = 12.5MK the middle of the 
temperature range. We estimate the volume using V = whL, 
with w the width, h the line-of-sight component and L the ex- 
tent perpendicular in the image to w. This is found using a filling 
factor of unity, and as this could be smaller, the values we cal- 
culate are upper limits to the thermal energy. Note that for the 
envelope region we have calculated the volume of the region mi- 
nus the volume of the plasmoid core. The energies are shown in 
Table [T] for the times and regions in Figure 12 We find that the 
thermal energy in this temperature range as the core erupts in- 
creases slightly (4.1 to 5.1 X 10^^ erg) whereas the energy in the 
stem structure greatly increases just after the impulsive eruption 
(1.2 to 9.1 X 10^^ erg). This is due to the widening of the stem 
whilst maintaining the same emission measure per pixel within 
the stem. This weakens by 12:18:02 UT when the emission mea- 
sure starts to drop and the stem narrows. If we again assume the 
structures are elongated, with line-of-sight h = lOOMm, we get 
a smaller variation between the thermal energy of the different 
structures and how they evolve with time, with energies of about 
^q30 gj.g ^Q^i^ cases there is an increase in the thermal energy 
of the stem which is consistent with the idea of continuous en- 
ergy release in this region. The envelope of hot material around 



the plasmoid core has a considerably higher thermal energy but 
that is because it is larger structure. 

We can also calculate the kinetic energy for the structures 
which clearly propagate outwards, the plasmoid envelope and 
core. The kinetic energy is taken as = nVm^v^ 12 where 
is the proton mass and the velocity is v = 660 kms"^ from Figure 
|9l These are again given in Table [T] for /z = w and h = lOOMm. 
Here we find that the thermal energy (in 11-14MK) is consis- 
tently slightly larger than the kinetic energy of the material in 
this temperature range. A higher thermal to kinetic energy was 
also found previously by'L andi et al.| ( |201 ). This similarity be- 
tween the thermal and kinetic energies is typical for CMEs from 
observations higher in the corona (e.g. |ForbeSj2 000). Our values 
for the ener gies are also s maller than those for a "moderately" 
large CME ( |Forbes||2000t (10^^ versus 10^^ erg) which is ex- 
pected as our energy is just from one component of the CME 
and in one temperature range. 

3.5.1. Coronal Loop Density 



In the later time intervals in Figure 12 (12:18:02 UT) there 
are distinctive coronal loops associated with the flaring region 
(LI) and the shrinkage below the eruption (L2). Again assum- 
ing w = hwe can calculate the loop density, taking them to be 
about 2Mm, though we do not estimate the energy as we cannot 
determine a reliable loop length. The large coronal loop (LI), 
associated with the core of the occulted flaring region, is visible 
from 12: 17:02 UT onwards, and its density increases as more hot 
material fiUs the loop (from 5.6x lO^cm"^ to 7.1 x 10^ cm"^). The 
loops below the ejected material (L2) have an emission measure 
that doubles between 12:17:02 UT and 12:18:02 UT, resulting in 
a density increase to 1.3 x 10 cm" from 8.9 x lO^cm'^ This 
increase in heated material maybe due to outflows from the re- 
connection region in the stem structure ( Savage et al. 2012) or 
the reconnection layer above the loops where Kelvin-Helmholtz 
waves were located ( [FouUon et al.|2Q12] ). 



4. Discussion & Conclusions 

Regularized inversion of multi-filter observations provides an es- 
sential tool to study the temporal and spatial evolution of multi- 
thermal plasma during eruptive events. The regularized maps of 
an erupting CME presented in the paper allow the basic plasma 
properties (integrated along the line-of-sight) to be inferred and 
calculation of the density and energetics of various parts of the 
eruption. We find that the leading edge of the eruption travels 
between 1 029-1 246 km s"^ and has the temperature range <2- 
3MK followed by slower 597-664 km s"^ but hotter plasma >3- 
4MK. The erupting core and stem which could be interpreted as 
a plasmoid with a current sheet appears sharply defined in 11- 
14MK. The width of the stem/current sheet we find (7-25Mm) is 
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far larger than that expected in the classical reconn ection model, 
and somewhat smaller than the predicted < I Mm CVrsnak et al.' 
|2QQ9j . This value is in fact smaller than those found higher in the 
corona (e.g. Ciaravella et a l. 2002; Lan di et al.|2012 ) suggesting 
continuous expansions as the CME lifts up. The stem/current 
sheet width we determine is however at a higher temperature 
range and made with a higher spatial resolution observations 
than these previous results. 

During the impulsive phase, we find that the density is about 
3 X 10^ cm"^, 6 X 10^ cm"^, 9 x 10^ cm"^ in the plasmoid core, 
stem and surrounding envelope of material. This gives thermal 
energy estimates of 5 x 10^^ erg, 1 x 10^^ erg and 2 x 10^^ erg. 
These are slightly larger than the values found for the kinetic en- 
ergy of the erupting envelope and core, which is consistent with 
previous CME energy estimations ( [Landi et"ar||2010| ). The ob- 
servations also show the increase of the thermal energy of the 
core (plasmoid) and of the stem (current sheet) as the CME is 
rising from 12:14:02 UT to 12:17:02 UT, suggesting ongoing 
energy release. The increase of energy is more noticeable for the 
stem with it more than doubling over 3 minutes. This is consis- 
tent with the predictions of recent reconnection models ( Murphy 
|et a l. 2010; Reeves et al.|2 010, e.g.) with the energy released in 
the current sheet pumped into the plasmoid. 

The reconstructed DEM maps show the multi- thermal nature 
of the corona, which is unsurprising given that we have a hot 
plasmoid erupting through the cooler surrounding corona, two 
distinct structures with diff'erent temperature characteristics. Our 
DEM maps show in greater detail the same basic physical picture 
that was found by [Cheng et al. (2012 ) using a diff'erent DEM re- 
construction method for a few regions in the CME. Namely that 
the bright leading front of the CME is composed of a corona- 
like temperature distribution, with the increased emission due to 
the accumulation of material, whereas the erupting core contains 
heated material (> 8MK). All DEM solutions however should 
be treated carefully given the uncertainties in the temperature 
resp onse functions an d the possibly broad temperature resolu- 
tion ( [Testa et al.|2012[ ). Used carefully, with these issues kept in 
mind, regularized DEM maps present a useful tool for explor- 
ing the dynamics of heating in the solar corona with SDO/AIA 
observations. 
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Appendix A: EM Movie of the Eruption 




Fig. A.l. One frame of the movie, available online, showing the 
EM maps, logn^h in units of cm"^, for the temperature ranges 
1-1.5 through to 14-19MK, starting at 12:12:02 UT running till 
12:19:52 UT. 



